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Abstract 

In present paper we propose seemingly new method for finding solu- 
tions of some types of nonlinear PDEs in closed form. The method is 
based on decomposition of nonlinear operators on sequence of operators 
of lower orders. It is shown that decomposition process can be done by 
iterative procedure(s), each step of which is reduced to solution of some 
auxiliary PDEs system(s) for one dependent variable. Moreover, we find 
on this way the explicit expression of the first-order PDE(s) for first inte- 
gral of decomposable initial PDE. Remarkably that this first-order PDE 
is linear if initial PDE is linear in its highest derivatives. 

The developed method is implemented in Maple procedure, which can 
really solve many of different order PDEs with different number of in- 
dependent variables. Examples of PDEs with calculated their general 
solutions demonstrate a potential of the method for automatic solving of 
nonlinear PDEs. 



1 Introduction 

Nonlinear partial differential equations (PDEs) play very important role in many 
fields of mathematics, physics, chemistry, and biology, and numerous applica- 
tions. If for nonlinear ordinary differential equations (ODEs) one can observe 
incontestable progress in their automatic solving, the situation for nonlinear 
PDEs seems as nearly hopeless one. 

Despite the fact that various methods for solving nonlinear PDEs have been 
developed in 19-20 centuries as the suitable groups of transformations, such as 
point or contact transformations, differential substitutions, and Backlund trans- 
formations etc., the most powerful method for explicit integration of second- 
order nonlinear PDEs in two independent variables remains the method of Dar- 
boux PQ-[1]. The original Darboux method (as already Darboux stated in pQ) 
is extendable in principle to equations of all orders in an arbitrary number of 
independent variables, even to systems of equations; however, in [T]-[2] and sub- 
sequent papers by many authors, the detailed calculations were performed only 
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for a single second-order equation with one dependent and two independent 
variables. 

The Darboux method was refined in recent years into more precise and effi- 
cient (although not completely algorithmic) form 5J-|8 and references therein. 
Nevertheless this approaches suffer from high complexity and necessitate to use 
some tricks. 

There were some partially successful attempts to extend modern variants of 
the Darboux method based on Laplace cascade method on higher-order PDEs 
and PDEs in the space of more than two independent variables [TU]-[T3] but 
they suffer from high complexity too. 

There is an original approach to the problem, based on the special type 
of local change of variables which leads to the order reduction of initial PDE, 
proposed in [14] . which is suitable for high dimensions problems but of very 
special class though. 

In present paper we propose seemingly new method for finding solutions 
of some types of nonlinear PDEs in closed form. The method is based on 
decomposition of nonlinear operators on sequence of operators of lower orders. It 
is shown that decomposition process can be done by iterative procedure(s), each 
step of which is reduced to solution of some auxiliary PDEs system(s) for one 
dependent variable. Moreover, we find on this way the explicit expression of the 
first-order PDE(s) for first integral of decomposable initial PDE. Remarkably 
that this first-order PDE is linear if initial PDE is linear in its highest derivatives. 

The developed method is implemented in Maple procedure, which can really 
solve many of different order PDEs with different number of independent vari- 
ables. Examples of PDEs with calculated their general solutions demonstrate a 
potential of the method for automatic solving of nonlinear PDEs. 



2 Bases of the method 

2.1 Decomposable PDEs 

The simplest second-order non-linear PDE for w = w(t,x) 

d 2 w 1 dw dw 
dtdx w dt dx 

can be easily transformed to the following decomposed form 



(1) 



izM~) = o, (2) 
at w ox 

from which we can without difficulty obtain the general solution to PDE |TJ) in 
two steps. First step gives us 

ldw = d\n(G(x)) 
w dx dx 



2 



where G{x) is an arbitrary function. And then, solving the equation (J3]) on the 
second step, we obtain 

w(t,x)=F(t)G(x), (4) 

where F(t) is one more arbitrary function. 

The main observations on analyzing the grounds of solvability of the PDE 
([I]) by the above method are that 

1. The PDE fTJ) is "decomposable", i.e., it can be represented as a composi- 
tion of successive differential operators of type (not necessarily linear) . It is 
clear that such type of decomposition can be done for some PDEs of any order 
and with any number of independent variables in the following manner 

Di(w) = ui , 
D 2 (ui) = u 2 , 

(5) 

-D„(u n -l) = 0, 

where x — (x\, . . . , x m ), w — w(x), Ui = Ui(x) and 

du du 
D i (u) = V i (x,u _,...,_). 

Assuming that Vi are arbitrary functions, and eliminating Ui by successive sub- 
stitutions in system |[SJ), we get a family of PDEs for w of nth order 

D n (D n - 1 (...D 1 (w)...))=0. (6) 

which are " decomposable" and in principle their solutions general or particular 
can be obtained by integration of split system ([5]). The PDE © is nonlinear if 
at least one of the operators Di is nonlinear. Not all PDEs admit such repre- 
sentation. And in positive cases such representation is not unique in general. 

Note that as a matter of fact Di need not be the first-order differential 
operators. So the composition procedure for nth order PDE, when n > 2 can 
be as follows 

D?(u)=0, (7) 
where ni, n 2 are integers and n\ + n 2 = n, w = w(x), u = u(x), and (k < j) 

DHu) V(x u — — —) 

OXi dX 1 ...dX r ^ n \k 1 +-+k m =k<j OXm 

The late representation allows us to carry out the PDE's decomposition or 
order reduction gradually bit by bit. 

We have to stress here that in general representations ([5]) and ([7]) may have 
different meaning. For example, some PDEs do not admit representation j5]) 
but permit the form ([7|) with both solvable DEs. 
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2. Each step of the solving process for decomposed PDE is faced with the 
necessity to solve differential equation Di{ui-\) — Ui (or D\ (itj_i) = Ui), so all 
such DEs must be solvable. Note that only first step D n (u n ^i) — is free from 
arbitrary functions. 

So one of the PDEs solving strategies may be as follows. First of all we 
try to decompose given PDE. In order to do so we have to solve corresponding 
auxiliary nonlinear PDE system for unknown functions Vi, it is sufficient to 
find a particular solution here. And, if it is successful, then, deciding between 
the variants, try to solve each arising DE from the chain ([5J. Main obstacle 
here, beginning at the second step is just mentioned necessity to solve DEs with 
arbitrary functions. There are sufficiently narrow circle of solvable (in sense of 
the general solutions) DEs with an arbitrary function as a parameter. 

Another (classification) approach can be based on the usage of only solvable 
DEs. That is, we can form a composition of successive solvable differential 
operators and as a result obtain a families of solvable PDEs. Such a way leads to 
extensive nontrivial families for different types of nonlinear PDEs which general 
solutions can be expressed in closed form. But on this way we encounter a 
difficulty to circumscribe such families integrally and are forced to consider 
particular subfamilies. Nevertheless it yields extensive field of PDEs for methods 
testing [15] . 



2.2 Decomposition algorithm for decomposable PDEs 

For nth order PDE, when n > 2 there are some slightly different approaches 
which are dictated by goals of the problem. If the goal is to decompose given 
nonlinear operator then we have to use the scheme ([7]) with n\ = 1, ri2 = n — 1. 
And conversely we have to use the scheme ([7} with m = n — 1, n% = 1 if the 
goal is to solve given PDE. The last procedure in some features resembles the 
well-known technics of reducing ODEs order, e.g., by first integral method. Of 
course, it is possible to use intermediate cases. 

All above cases can be treated by the same way as we consider below but 
each of them leads to auxiliary PDEs systems of different order, viz n^ + l, with 
corresponding calculation complexity. 

In sequel we will consider for shortness only the case with ri\ = n— 1, n<i = 1, 
as more practical for PDEs solving. 

Let us consider the decomposition of type ((Jj with D™~ l (w) as a solution 
of the following equation with respect of u = u(x) 

J{n,x,w,— r — — „_ 1 = 8) 

and 

D 2 (u) = V(x,u,-^,...,^). (9) 
0x1 ox m 

If substitute u = D"^ 1 (w) into © we obtain decomposable n-th order PDE 

V(x,U ,U Xl ,...,U Xm ) = 0, (10) 
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where (we use below the following notation w = Wo and — k f w k — Wk 

dx x . ..dx^ 

~\7l—l I 



Dr i N=f/ , (ii) 

dJ 



dJ 
On 



= U Xi (i = l,...,m), (12) 



where k* = fcj + 1 if j = i and k* — kj otherwise, and it is supposed that 
differentiations in sum are carried out on all indexed W : s which are involved in 
J. 

Here we can introduce Uq and U X1 , . . . , U Xm as new independent variables 
if express m variables from the set {Wk lt ....k m } with k\ + • • • + k m = n using 
linear system ([12")) . 

Assuming that given PDE of order n 

dw d k w d n w 

is decomposable, we receive, that after substitution of the new variables, left- 
hand side of given PDE must turn into (|10[) with some V. 

Left-hand side of given PDE expressed in new variables is the first-order 
differential expression with respect to 

J(U ,x, W , Wi t o,...,o, ■ ■ ■ , Wfe 1 ,...,fc m |fe 1 +...+fc m =fe<n-i, ■ • ■ , Wo,o,...,n-i) 

and must not depend on all indexed W's, that is derivatives of F expressed in 
new variables with respect to all indexed W's are equal to zero. Sequence of 
such derivatives of F equated to zero form a second-order PDE system for J. 
So a solution (particular as well) the PDE system gives possible expression of 
differential operator D r l~ 1 {w) through (jHJ) and differential operator D 2 (u) by 
substituting the solution of J into left-hand side of given PDE expressed in new 
variables. 

Of course, there are problems where a operator decomposition is required 
only. But in most cases obtained decomposition is intended for finding solutions 
for given PDE. If in obtained decomposition the corresponding PDE D2(u) = 
is solvable, then substitution of obtained u into J expressed in original variables 
gives us a first integral (see its definition in the next subsection) of given PDE. 
It is easy to see that for decomposable PDEs the first integral is a differential 
equation, so we can try to solve it or to find a first integral for this new DE (or 
decompose it) by the scheme described above until we come to the first-order 
DE. 

Remarkably that in the approach under consideration the finding of first 
integrals can be done more directly and effectively. 



2.3 Differential equation for first integral of decomposable 
PDEs 

The first integral / of the PDE is an expression, involving one arbitrary func- 
tion, which is equivalent in some sense to the given PDE. The first integral 



5 



vanishes on the set of solutions of given PDE. And (in accordance with [1]) all 
differential consequences of the equation 1 = coincide with respective differen- 
tial consequences of given PDE (e.g., elimination of the arbitrary function leads 
to the given PDE). 

Our goal here is to find PDE for first integral of a decomposable PDE. To 
do so we first of all have to take into account that u{x) is the solution of the 
corresponding PDE 

, t du du . 

^ (x ' U 'fe?'--'c^ ) = ' 

so u(x) depends only on x but in no way on indexed W l s. Secondly, the depen- 
dent variable in this case, namely 

J{u(x),X, Wq, Wl,0,...,0, • ■ ■ , Wk u ...,k m \k 1 +-+k m =k<n-l, • ■ ■ , Wo,0,...,n-l) 

of given PDE l|13p expressed in new variables do not to depend on U X1 , . . . , U Xm 
and is a first integral of given PDE. 

If now consider u(x) as an unknown function, we can denote the first integral 

as 

I(x, W , Wl,0,...,0) ■ • • j Wf, u ...,k m \ki+-+k m =k<n-U ■ ■ ■ i Wo,0,...,ra-l) = 
J(u(x), X, Wq, Wl,o,...,0, ■ ■ ■ , Wk u ...,k m |fci+-+fc m =fc<n-l) • • • ) Wo,0,...,n-l) 

and instead of (fT2]) in the form 
dJ_ \ - 8J 

<)■>■, ' ^nw,, : 

we arrive to the following system 

e+Esb^x^ ('-'.-■■»)• <«> 

If express m variables from the set {Wkj, fc m } with k\ + • • • + k m = n (at 
least one of which is actual for given PDE - note that there are some variants 
here as a rule, so we can obtain some consistent PDEs on this stage) using linear 
system (|14|) and substitute them into given PDE (fT3|) we receive a first-order 
(even linear if PDE (|13p is linear in its highest derivatives) PDE with respect to 
first integral /. And it remains only to solve this PDE(s) to find a first integral 
of given PDE. 

Note, given PDE is decomposable iff exists a solution of such first-order 
PDE(s). 

3 Examples 

To facilitate necessary calculations in the process of finding first integrals I 
have implemented above described method in prototype of Maple procedure 
reduce_PDE_order (see Appendix). The input data of the procedure are given 



-U, 



dJ_ 

* du 



(i = \,...,m) 
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PDE of any order and dependent variable of the PDE with any number of 
independent variables. The procedure tries to find first integral(s) of the input 
linear or nonlinear PDE. 

The Maple built-in procedure pdsolve is used inside my procedure to solve 
the first- order PDE for first integral. As different Maple versions have different 
PDE solving abilities so the output depends on Maple version. In the following 
examples I refer to Maple 11. 

The procedure reduce_PDE_order is able to find first integrals for many 
known and unknown linear and nonlinear PDEs. Here we give examples of 
PDEs for which it is possible to find finally their general solutions. More exam- 
ples one can find in collection of solvable nonlinear PDEs [15 . 



3.1 Second-order PDE with two independent variables 

For PDE (w = w(t,x)) 



d 2 w a ( dw\ 
dtdx w \ dx ) 
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1 dw c \ dw 

w dt w I dx 2aw dt 



c dw be 
— kw— - 



2a 4aw 



(15) 



with o / and 4ak — b 2 ^ the procedure reduce_PDE_order outputs the 
following first integral 



I = Fl 



, ( 'c+2a^ 

t^jAak - b 2 - 2 arctan dx 

\ icv 4ak — b 2 



x, ■ 



V4afe - b 2 



with arbitrary function Fl. 

The ODE I = can be solved and one obtains (after some hand simplifica- 
tions and edition) the following general solution to (|15[) 



w(t, x) 



2a 



exp(Wb 2 - 4ak)F(x)(b + ^Jb 2 - Aak) - ^Jb 2 - Aak + b 
1 + cxp(Wb 2 - 4ak)F(x) 



dx 



dx+ 



G(t)} exp 



1 

2a 



exp(tV& 2 -4ak)F{x){b + Vb 2 - Aak) - Vb 2 - 4ak + b 



1 + cxp(Wb 2 - 4ak)F(x) 
where F(x) and G(t) are arbitrary functions. 



dx 
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3.2 Second-order PDE with four independent variables 

For PDE 



A 1 



d 2 w 



A, 



d 2 w 
' 8x28x4 



. 8 2 w „ 8w 
A 3 f . » +Co + ^i» — 
8x30x4 0x4 



8x18x4 

„ , . dw . 8w dw „ „ . 
C7i(Ai— + A 2 — + A 3 — + Biw + 5 ) + 

OTl 8X2 8X3 

CM^i — + A 2 — + A 37 — + B lW + Bo 2 = , 
8x\ 8x2 8x3 



(16) 



where w — w(xi,x%, X3, X4) and Ai,Bi,Ci are constants, the procedure re- 
duce_PDE_order outputs the following first integral 



I = Fl 



2 arctan 



X 1: X 2 ,X3,X 4 + 



2C7 2 (A 1 ^- + ^ 2 ^ + ^ 3 ^ + B iM , + Bq) + C 1 
\J 4C0C2 — C\ 
y/4C C 2 - cf 



with arbitrary function Fl. 

The PDE 1 = can be solved and one obtains the following general solution 

to iia 



tu(£i,a;2,ar3,X4) = 

cxp(-^i) ^ exp(^)(2i? C 2 +d + t<m[~x^4C C 2 - C? 



2Ai 



G(£, (A 2 £ + Aix 2 - A 2 a;i), (A 3 £ + ^3 - A 3Xl ))]J 4C C 2 - C 2 )d^ 



+ exp(- 



At 



-)F[(AiX2 - A 2 Xi), (A1X3 - A3X1), X4] , 



where F(ti,t 2 ,t3) and G(ti,t 2 ,ts) are arbitrary functions, c is arbitrary con- 
stant. 



3.3 Third order PDE with two independent variables 

For PDE (w = w(t,x)) 



8 3 i 



2w 



d 2 w dw n 8w ( 8w\ 



+ 2— — 

dtdx 2 dtdx 8x dt \dx J dt 8x 2 

the procedure reduce_PDE_order outputs the following first integrals 



dw 8 2 



aw 6 = 



(17) 



h = Fl 



1 , 8 2 w dw dw 



1 



w 2 dtdx dt dx 



axw 



ax 2 w 2 + 2w( 



dw 
~dt 



d 2 w 
'dtdx' 



2x 



dw dw 
dt dx 
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and 



h = Fl 



dx 2 



atw 



dw\' 
dx J 



with arbitrary function Fl. 

We can form some PDEs from Ii and to solve them we can repeat the process 
of order reduction with the procedure reduce-PDE_order. The ODE 1% = can 
be solved directly and one obtains in any way the following general solution to 

m 

a(t , x) _ mmr (ff _ xm +x J G(x)dx _ j xG(x)ix ) , 

where F(t), H(t) and G(x) are arbitrary functions. 

3.4 Fourth order PDE with two independent variables 

For PDE (w = w(t,x)) 



, d 4 w „ 9 / d 3 w dw d 3 w dw 

W A 2w I 

dt 2 dx 2 \ dt 2 dx dx dtdx 2 dt 



2 w 



d 2 w 
dtdx 



} dw dw 
J dt dx 



dx 2 



dw 
dx 



dw 

at 



dt 2 



dx 2 



dx 



= 



the procedure reduce_PDE_order outputs the following first integrals 



(18) 



h=Fl{t, 



d :i 



w 



w 



dt 2 dx 



2w 



dw d 2 w 
dt dtdx 



d 2 w dw n dw ( dw 



dt 2 dx 



dx V dt 



d\ 



dt 2 



dt 2 dx 



dw d 2 w 



dt dtdx dt 2 dx 



dw 
~dt 



dw ( dw\ 

w - 2x dx~ {-m) 



and 



h = Fl(x, ■ 



2 d 3 w ^ 9w d 2 w dw d 2 w ^dw / dw 
dtdx 2 dx dtdx dt dx 2 dt \ dx 



1 



d 2 w 



t- 



d 3 w 



dx 2 dtdx 2 J 
with arbitrary function Fl. 



dw d 2 w d 2 w dw 



dx dtdx dx 2 dt 



( dw 
V dx 



w-2t 



dw ( dw 



dt \ dx 
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The wealth of first integrals here allows us to operate with them in many 
different ways. Apart from aforesaid subsequent order reduction we can, for 
example, from 



1 



, d 3 w ^ 9w d 2 w dw d 2 w ^dw ( dw 
dtdx 2 dx dtdx dt dx 2 dt \ dx 



Fix) 



and 



1 J d 2 w ^ d 3 w \ 2 
dx 2 dtdx 2 J 



^dw d 2 w d 2 w dw\ ^ ( dw 
dx dtdx dx 2 dt J \ dx 



where F(x) and G(x) are arbitrary functions, algebraically eliminate mixed 
derivative and obtain the following ODE 



dw 



dx 2 



+[tF(x)-G(x)]w 2 = 0, 



which gives the general solution to ([18 



j(t, x) = H(t) exp 



t / xF(x) dx - tx / F(x) dx+ 



x J G{x) dx- J xG(x) dx + xK(t) 
where F(x), H(t), G(x) and K(t) are arbitrary functions. 



4 Conclusion 

The method have considered above is efficient enough for solving decomposable 
PDEs of relatively high order with many independent variables. The main 
limitation here is concerned with abilities to solve corresponding auxiliary first- 
order PDEs for first integrals. 

An adaptability of the method to PDEs which are not decomposable but 
which general solutions can be expressed in closed form remains unsolved yet. 
But it can be shown on examples that there are some ways to extend the method 
for some types of such PDEs. These approaches deserve further thorough study 
in another publication. 
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5 Appendix. 

Maple procedure reduce_PDE_order 

rcducc_PDE_ordcr : =proc (pde,unk) 

local B,W,N,NN,ARG,acargs,i,M,pde0,DN,IND,IND2,IND3,IND4,ARGS,SUB,SUB0, 
ZO,Bargs,EQS,XXX,WW,BB,PP,pdeI,IV,s,AN; 

option 'Copyright (c) 2006-2007 by Yuri N. Kosovtsov. All rights reserved.'; 
N:=PDETools[difforder](op(l,[selectremove(has,indets(pde, function), unk)])); 
NN:=op(l,[selectremove(has,op(l,[selectremove(has,indets(pde, function), unk)]),diff)]); 
ARG:=[op(unk)]; 
acargs: ={}; 

for i from 1 to nops(ARG) do 

if PDETools[diff order] (NN,op(i,ARG))=0 then else acargs:=acargs union {op(i,ARG)} 
fi; od; 

acargs:=convert(acargs,list); 

M:=op(0,unk)(op(acargs)); 

if type(pde,equation)=true then 

pdeO : =lhs(subs(unk=M,pde) )-rhs (subs (unk=M,pde) ) else pdeO : =subs (unk=M,pde) 
fi; 

DN:=[seq(seq(i,i=l..nops(acargs)),j=l..N)]; 
IND:=seq(op(combinat[choose](DN,i)),i=l..N); 
IND2:=seq(op(combinat [choose] (DN,i)),i= 1 . .N-2) ; 
IND3 : =op (combinat [choose] (DN,N- 1 ) ) ; 
IND4:=op(combinat[choose](DN,N)); 

ARGS:=op(unk),M,seq(convert(D[op(op(i,[IND2]))](op(0,unk)) 
(op (acargs) ) ,diff ) ,i= 1 . .nops ( [IND2] ) ) ; 

SUB:={M=W[0],seq(convcrt(D[op(op(i,[IND]))](op(0,unk)) 

(op(acargs)),diff)=W[op(op(i,[IND]))],i=l..nops([IND]))} ; 

SUB0:={W[0]=op(0,unk)(op(ARG)), 

seq(W[op(op(i,[IND]))]=subs(M=op(0,unk)(op(ARG)), 

convert(D[op(op(i,[IND]))](op(0,unk))(op(acargs)),difF)),i=l..nops([IND]))}; 

Z0:=B(ARGS,seq(convert(D[op(op(i,[IND3]))](op(0,unk))(op(acargs)),diff), 

i=l..nops([IND3]))); 

Bargs: =op (indcts (subs (SUB , Z0) ,name) ) ; 

EQS:=convert(subs(SUB,{seq(diff(ZO,op(i, acargs) )=0,i=l..nops(acargs))}),diff); 

XXX:={seq(W[op(op(i,[IND4]))],i=l..nops([IND4]))}; 

WW:=select(type,indets(subs(SUB,pdcO)), 'name') intersect 

{seq(W[op(op(i,[IND4]))],i=l..nops([IND4]))} ; 

BB : =sclect (has, combinat [choose] (XXX, nops (acargs) ) , WW) ; 

PP:={}; 

pdeI:={seq({subs(subs(solve(EQS,op(i,BB)),subs(SUB,pdeO)))},i=l..nops(BB))}; 
IV: ={seq( W [op (op (i, [IND4] ) )] ,i=l . .nops( [IND4] ) ) } ; 
for s from 1 to nops(pdcl) do 
try 

AN:=pdsolve(op(s,pdeI),{B},ivars=IV); 
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for i from 1 to nops(AN) do 
if op(0,lhs(op(i,AN)))=B then 
PP:=PP union {rhs(op(i,AN))} 

fi; 

od; 
catch: 
end try; 
od; 

PP:=subs(SUBO,PP); 
RETURN(PP); 
end proc: 



Calling Sequence: reduce_PDE_order (PDE, f{x)); 

PDE - partial differential equation; 

f(x) - indeterminate function with its arguments. 
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